The Shear Viscosity to Entropy Density Ratio 
of Trapped Fermions in the Unitarity Limit 
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We extract the shear viscosity to entropy density ratio r//s of cold fermionic atoms in the unitarity 
limit from experimental data on the damping of collective excitations. We find that near the critical 
temperature r\j s is roughly equal to 1/2 in units of h/ks- With the possible exception of the quark 
gluon plasma, this value is closer to the conjectured lower bound l/(47r) than any other known 
liquid. 



I. INTRODUCTION 

Strongly correlated quantum systems play an impor- 
tant role in many different areas of physics. Systems of 
interest in atomic, condensed matter, and nuclear physics 
span many orders of magnitude in energy scale but ex- 
hibit a number of universal properties. Recently there 
has been renewed interest in the transport properties of 
strongly correlated systems. Experiments at the Rela- 
tivistic Heavy Ion Collider (RHIC) indicate that at tem- 
peratures close to the critical temperature T c the quark 
gluon plasma is strongly interacting. The strongly inter- 
acting quark gluon plasma (sQGP) is characterized by a 
very small shear viscosity to entropy ratio, a small heavy 
quark diffusion coefficient, and a large opacity for high 
energy jets 0, 1 1 II. 

From a theoretical point of view not much is known 
about transport coefficients of strongly correlated sys- 
tems. If the interaction is weak then the mean free path 
and the shear viscosity are large. As the strength of the 
interaction increases the mean free path and the viscosity 
drop but there are good reasons to believe that the shear 
viscosity always remains finite. Kovtun et al. conjectured 
that there is a universal lower bound r\j s > %/ (inks) Q. 
Here, r\ is the shear viscosity, s is the entropy density, h 
is Planck's constant and ks is the Boltzmann constant. 
The bound is saturated in the case of strongly coupled 
gauge theories that have a dual description in terms of a 
gravitational theory. 

In this work we test the viscosity bound conjecture 
by extracting ry/s from experimental data on the damp- 
ing of collective oscillations of a cold atomic gas near a 
Feshbach resonance. Cold atomic gases provide an ideal 
system to test the conjecture because both the temper- 
ature and the interaction can be continuously adjusted. 
Also, because of universality, atoms in the unitarity limit 
are equivalent to other Fermi liquids with a large scatter- 
ing length like dilute neutron matter. Collective modes 
in the atomic system have been studied in a number of 
experiments 0, H, 0, EH- In the weak coupling regime 
the frequency and damping constant of collective modes 
can be understood in terms of the Boltzmann equation. 
In the unitarity limit the frequency of collective modes 
is well described by ideal hydrodynamics [III E2, EH ■ In 
the present work we include viscous corrections and use 



experimental and quantum Monte Carlo data on the ther- 
modynamics to extract rj/s. Previous studies of damping 
near the unitarity limit can be found in [3, EH, EH • 



II. EULER HYDRODYNAMICS 

We shall assume that the system is approximately de- 
scribed by ideal (Eulerian) fluid dynamics. The equation 
of continuity and of momentum conservation are given 
by 



dn 
~dt 



+ V • (nv) = 0, 



dv 



mn— + ran ( v ■ V ) v = —VP — riVV, 
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where n is the number density, m is the mass of the 
atoms, v is the fluid velocity, P is the pressure and V is 
the external potential. The trapping potential is approx- 
imately harmonic 



2 2 
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In the unitarity limit the equation of state at zero tem- 
perature is a simple polytrope P ~ n 7+1 with 7 = 2/3. 
At finite temperature the equation of state is more com- 
plicated, but universality implies that the isentropic com- 
pressibility is unaffected, 



dP 

dn 



(7 + 1) 



P 
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The equilibrium distribution no can be determined from 
the hydrostatic equation VPq = — rioVV. At T = 



n (f) =n (0) 1 
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where /i is the chemical potential. In the unitarity limit 
the chemical potential is related to the Fermi energy as 
fi = £,Ep, where £ is a universal parameter (£ ~ 0.44 
according to the quantum Monte Carlo calculation [TtJ )• 
The central density and the total number of particles are 



N = 



1 



3£ 3 /2 
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FIG. 1: Temperature dependence of the damping rate for the 
radial breathing mode of a trapped 6 Li gas near the Feshbach 
resonance at 840 G, from Kinast et al. [jj. Here r is the 
damping time and ui± is the radial trap frequency. We have 
used the calibration curve in pl| to convert the f scale of [H 
to T/T F . 

where to — {uj\uj2Lo^) 1 ^ ■ Consider small oscillations 
n = riQ + Sn. From the linearized continuity and Eu- 
ler equation we get [Tl| 

8 2 v 



dt 2 



-1 



(v-u) (vv) -V(V-W) 
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where we have dropped terms of the form Vi'Vjif that 
involve higher derivatives of the velocity. Using a scaling 
ansatz vi — aiXi exp(zu;i) (no sum over i) we get 



(2a, 



oj 2 ) a. 
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This is a simple linear equation of the form Ma = 0. Non- 
trivial solutions correspond to det(M) = 0. In the case of 
a trapping potential with axial symmetry, qjj — ll>2 = <-^q, 
lu 3 = Xu; , we get uj 2 = 2cv 2 and [ill, El, Ell 



to 



^ 7 + l + I+^A 2 
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In the unitarity limit (7 = 2/3) and for a very asymmetric 
trap, A — » 0, the eigenfrequencies are uj 2 — 2w 2 and 
lo 2 = (10/3)w 2 . The mode u 2 = (10/3)w 2 is a radial 
breathing mode with a = (a, a, 0) and the mode uj 2 = 2cJq 
corresponds to a radial quadrupole a = (a, —a, 0). 



III. VISCOUS CORRECTIONS 



The energy dissipated due to viscous effects is 
E = ~ I d 3 x7i(x) (diV-j+djVi-^SijdkVk) (10) 



d 3 x({ X ) (diVif 



where i] is the shear viscosity and £ is the bulk viscos- 
ity. In the unitarity limit the system is scale invariant 
and the bulk viscosity in the normal phase vanishes. The 
situation is more complicated in the superfhiid phase. 
In this case the normal and superfluid components can 
flow independently and in addition to the shear viscos- 
ity there are three bulk viscosities Q. Two of the three 
bulk viscosities, £i and £2, can be shown to vanish as a 
consequence of scale invariance (i~9j . The third bulk vis- 
cosity, £3, only contributes to dissipation if the normal 
components is moving relative to the superfluid [20| . In 
the following we shall consider modes in which the two 
components move in phase and ignore bulk viscosity. We 
will also assume that the viscosity only depends on x 
through the local density and temperature. This is valid 
if the density and temperature are varying slowly, and is 
consistent with the local density approximation for the 
density profile. We note that in this approximation the 
equation of state and the transport coefficients reflect the 
conformal invariance of the microscopic dynamics, even 
though scale invariance is broken by the external poten- 
tial. 

For the radial scaling flows given in equ. © we have 



E 



(al + a 2 , - a x a y ) d 3 xr](x) 7 (11) 



where E is a time average. The damping rate is deter- 
mined by the ratio of the energy dissipated to the total 
energy of the collective mode. The kinetic energy is 



rn 
T 



d 3 x n(x)v 2 



mN 



'a 2 x + a 2 y )(x 2 ). (12) 



At T = we find (a; 2 ) = R 2 j_/8, where R±_ is the trans- 
verse size of the cloud. At non-zero temperature we can 
use the Virial theorem [2l[ to relate (a; 2 ) to the total en- 
ergy of the equilibrium state, (x 2 ) / (x 2 )t=q = E/Et=o- 
The damping rate is [TU, HI] 



1 E _ 2 a 2 + a 2 - a x a y J d 3 x r)(x) 
2E ~ 3 a 2 +a 2 mN{x 2 ) 



(13) 



Note that the second factor on the RHS is 1/2 for the ra- 
dial breathing mode and 3/2 for the radial quadrupole 
mode. If this dependence could be demonstrated ex- 
perimentally, it would confirm that the damping is in- 
deed dominated by shear stress. Another possibility is 
to compare the breathing mode with a scissors mode. 
The scissors mode is characterized by the velocity field 
V = a\7(xy). The frequency is lu 2 = u> 2 
second factor in equ. (fl"3"| is 6. 

For the unitary Fermi gas the ratio of the shear vis- 
cosity to the entropy density is given by universal func- 
tion that depends only on the ratio T/Tp, rj(fi,T) — 
a(T/Tp)s([i,T). Here, the Fermi temperature is given 
by T F = (37r 2 n) 2 / 3 /(2m). In the local density approxi- 
mation this implies that rj(x) = a(T/Tp(x))s(x), where 
Tp{x) is the local Fermi temperature. We shall assume 
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FIG. 2: Viscosity to entropy density ratio of a cold atomic 
gas in the unitarity limit. This plot is based on the damping 
data published in [8j and the thermodynamic data in [ia. |23|. 
The dashed line shows the conjectured viscosity bound rj/s — 
l/(47r). The shaded band is a systematic error estimate based 
on the contribution to E from atoms outside a surface at 
optical depth one. 



that a is a smooth function and replace Tp(x) by its 
value at the center of the trap. This approximation can 
be checked a posteriori. We note that since the flow 
profile has a simple scaling form the damping rate is pro- 
portional to the volume integral of the shear viscosity. If 
t] ~ s then the damping rate scales with the total entropy. 
The kinetic energy, on the other hand, is proportional to 
the number of particles. The shear viscosity to entropy 
density ratio extracted from the radial breathing mode 
is 



n = ±e /2 (3N) 1/3 

s 4 
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E 
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(14) 



where T/lo±_ = 1/(tu>±) is the dimensionless damping 
rate. 



IV. BREAKDOWN OF HYDRODYNAMICS 

Near the edge of the cloud the mean free path is com- 
parable to the characteristic length scale over which the 
velocity field varies and the hydrodynamic description 
breaks down. In this regime the density is low and the 
mean free path can be estimated using kinetic theory. In 
the limit n <C (mT) 3//2 the mean free path is given by 



Imfp = { n(J Y 



4tt 1 
3 mksT 



(15) 



Following [22| we can define a surface r Q (8) by the con- 
dition that a particle incident from the exterior of the 
cloud has a chance of no larger than 1/e of colliding with 
another particle 

1=1 (16) 



Here, 9 denotes the angle between ds and the z-axis. 
Ideally, the hydrodynamic description inside the surface 
ro{0) should be coupled to kinetic theory outside this 
surface. In this work we shall be less ambitious and use 
equ. (|16[) in order to estimate the systematic uncertainty 
in our determination of rj/s. For this purpose we compute 
the contribution to E in equ. (|13jl that comes from atoms 
outside ro(0) and treat it as a systematic error in the 
damping rate. 



V. THERMODYNAMICS 

Experimental results for T/ui are shown in Fig. fl}. 
In order to extract rj/s we also need information on the 
entropy per particle. In the unitarity limit there are only 
two energy scales in the problem, the chemical potential // 
and the thermal energy fc^T (in the following we will set 
ks to one). The associated density scales are nf(ii) and 
A^ 3 , where rif(fx) = (2m/i) 3 / 2 /(37r 2 ) is the density of a 
free Fermi gas and A T = {ZirmT)- 1 / 2 is the thermal wave 
length. All thermodynamic quantities can be expressed 
as suitable powers of either n/(/z) or A?r times a function 
of the dimensionless quantity y = T/fi. For example, we 
can write the pressure as [24j 



P(ji,T) = |/in/0*)S(») = ^TX^Wiy- 1 ) , (17) 

where the first form is more useful at small T and the sec- 
ond at high T. Using standard thermodynamic identities 
one can show that 



n(/i,T) = n f (p)f{y), 
2 

s(p,,T) = -n f (n)g'{y), 



(18) 
(19) 



where n is the density, s is the entropy density, and 
T( y ) = Q( y ) - 2yG'(y)/5. At T = the function 
9(y) is related to the parameter £ introduced above, 
g(0) = ,e- 3 / 2 . The functions G(y), T(y) refer to the 
bulk system. Trapped systems can be described using 
the local density approximation. Experiments typically 
involve 10 5 — 10 6 atoms and the local density approxima- 
tion is very accurate. The density of the trapped system 



n (x) =n(fi- V(x),T) , 



(20) 



r (9) 



Imfp ( r ) 



where V(x) is the trapping potential. Similar relations 
hold for the energy and entropy density. For T = 
equ. (p?TJ)) reduces to equ. ©■ 

The function G(y) can be extracted from quantum 
Monte Carlo data or from calorimetric experiments with 
trapped fermions. A number of Monte Carlo calculations 
have appeared over the last couple of years [lH, [H, [13] , 
but there are still significant disagreements between the 
results. Burovski et al. find T c /Tp = 0.15 2(7) and a 
critical entropy per particle S/N = 0.16(2) [20|. Bul- 
gac et al. quote T C /T F = 0.23(2) and S/N ~ 1.1 
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Using the local density approximation the results of Bul- 
gac et al. correspond to a critical entropy per particle of 
S/N ~ 2.1 for the trapped system. The reason that the 
entropy per particle is larger in a trapped system is that 
the density near the edges is smaller, and therefore the 
dimensionless temperature mT/n 2 ^ 3 larger. 

The Duke group has performed a series of calorimetric 
measurements [18|,|23(. Kinast et al. [18| provide a simple 
parametrization of the energy of the trapped system as 
a function of t = T/Tp. The result is 



E 



97.3 1 373 
4.98 1 1 - 43 



t < t c 
t > t c 



(21) 



with E Q = 0.53-Ef and T C /T F = 0.29(2). Luo et al. j2J 
give a similar parametrization of the entropy, 



S/N 



4.6 (e 
4.0 (e 



„ \0.61 
„ \0.45 

eoj 



e < e c 
e > e n 



(22) 



where e = E/Ep. The critical entropy per particle is 
S/N — 2.7, roughly compatible with the Monte Carlo 
results of Bulgac et al., but significantly larger than the 
results of Burovski et al. 



VI. RESULTS AND CONCLUSIONS 

In Fig. [2] we show rj/s extracted using both the damp- 
ing data and the calorimetry from the Duke group. We 
observe that rj/s is small (~ 0.5) near T c and slowly 
grows with temperature for T > T c . The value of rj/s 
near T c is about six times larger than the conjectured 
viscosity bound and consistent with the picture of a very 
strongly correlated liquid. Indeed, the extracted value 
•q/s ~ 0.5 is smaller than the previously known minimum 
for all other liquids, rj/s ~ 0.7 for liquid Helium near the 
lambda point [6j. Even smaller values of rj/s have been 
reported for the quark gluon plasma produced at RHIC, 
but the uncertainties remain large @, [28| . We note that 
the systematic uncertainty near the minimum is small, 
but that it increases as a function of temperature. 

There are a number of issues that need to be addressed 
in more detail. We argued that it is important to estab- 
lish that shear viscosity is indeed the dominant damping 
mechanism. This can be done either by studying the de- 
pendence of the damping time on the type of collective 



mode, or by studying the dependence on system size. Ki- 
nast et al. collected some data on system size dependence 
below T c and find that r/wj_ is roughly independent of 
the number of particles This is not consistent with 
the scaling in equ. (|14[) . Data on system size dependence 
can also be used to determine at what point hydrody- 
namics is breaking down. 

It is also important to understand the damping mecha- 
nism below T c in more detail. The data show a very sim- 
ple linear behavior in the variable f ~ (T/T F ) 2 / 3 . The 
natural framework for understanding the damping mech- 
anism in the regime below T c is superfluid (two-fluid) 
hydrodynamics [20j. There are several sound modes in 
superfluid hydrodynamics. First sound is an excitation 
in which the superfluid and normal components move to- 
gether, whereas second sound corresponds to an oscilla- 
tion of the superfluid component against the normal one. 
The damping of first sound is governed by the shear stress 
of the normal component. It is likely that the collective 
excitations that have been observed experimentally are 
ordinary (first) sound modes, but it is not obvious why 
the damping constant is linear in T. We should note that 
a linear behavior was observed in trapped Bose gases [2!| , 
where it was attributed to Landau damping by normal 
excitations [30(. 

Finally, there are a number of technical aspects of our 
analysis that should be improved. We have assumed that 
the quantity a = rj/s is only weakly temperature de- 
pendent. Near the the minimum of rj/s this is a good 
approximation, but at higher temperature the uncer- 
tainty inherent in this approximation grows, as does the 
uncertainty related to the breakdown of hydrodynam- 
ics near the the surface of the cloud. In this regime a 
Boltzmann description should be used. In the present 
work we have neglected dissipation due to temperature 
gradients and thermal conductivity. This is expected 
to be a good approximation for scaling flows because 
oscillations in density are proportional to the equilib- 
rium density Sn(x) ~ no(x). For isentropic oscillations 
ST ~ (Sn/n)T and to leading order no temperature gra- 
dients are present. 
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